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Abstract 

In this paper we study the optimal m-states switching problem in finite horizon as well 
as infinite horizon with risk of default. We allow the switching cost functionals and cost 
of default to be of polynomial growth and arbitrary. We show uniqueness of a solution for 
a system of m variational partial differential inequalities with inter-connected obstacles. 
This system is the deterministic version of the Verification Theorem of the Markovian 
optimal m-states switching problem with risk of default. This problem is connected with 
the valuation of a power plant in the energy market. 
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1 Introduction 

In this work we are concerned with the following systems of m variational partial differential 
inequalities with inter-connected obstacles: 

min{t>j(t, x) — I max{—gij{t,x) + Vj{t,x)) V {—Fi{t,x)) 

-dtVi{t, x) - Avi{t, x) - ^i{t, x)} = 0, V(t, x) e [0, T]x JR'', iel= {1, m}, 
^ Vi(T,x) = 0. 

(1.1) 
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mm{vi{x) - max(-fc(z) + Vj^x)) V (-Fi(x)) ,rvi{x) - Avi{x) - i{ji{x)} = 0, 
\lxeM^, iel= {1, ...,m}, 




(1.2) 



where gij, ipt and Fi are continuous functions, A an infinitesimal generator associated with 
a diffusion process and finally X~* := {1, i — 1, z + 1, m}. 

These system is the deterministic version of the Verification Theorem of the optimal multi- 
modes switching problem in finite horizon and infinite horizon with risk of default. These 
problems, of real option type, can be introduced with the help of the following example: 

Assume we have a power station/plant which produces electricity and which has several 
modes of production, e.g., the lower, the middle and the intensive modes. The price of electricity 
in the market, given by an adapted stochastic process {Xt)t<T, fluctuates in reaction to many 
factors such as demand level, weather conditions, unexpected outages and so on. On the other 
hand, electricity is non-storable, once produced it should be almost immediately consumed. 
Therefore, as a consequence, the station produces electricity in its instantaneous most profitable 
mode known that when the plant is in mode « G X, the yield per unit time is given by means 
of i/ji, switching the plant from the mode i to the mode j is not free and generates expenditures 
given by g^j and, on the other hand, cost of default (definitely stop the production) in mode 
i G X, is given by Fj and possibly by other factors in the energy market. So the manager of the 
power plant faces two main issues: 

{i) when should he decide to switch the production from its current mode to another one? 

(ii) to which mode the production has to be switched when the decision of switching is made? 

In other words she faces the issue of finding the optimal strategy of management of the plant. 
This issue is in relation with the price of the power plant in the energy market. 

Optimal switching problems for stochastic systems were studied by several authors (see e.g. 
PEliaiaEl[IIl[I3[ISl[I3IISl[IZl[ia[221EniE3]and references therein) . The motivations are 
mainly related to decision making in the economic sphere. Several variants of the problem we deal 
with here, including finite and infinite horizons, have been considered during the recent years. 
In order to tackle those problems, authors use mainly two approaches. Either a probabilistic one 
[T2| [T3t \T6\ [T7t 122] or an approach which uses partial differential inequalities (PDIs for short) 



The PDIs approach turns out to study and to solve, in some sense, the system of m PDIs 
with inter-connected obstacles ( II. ip for finite horizon and (11. 2p for infinite horizon. Recently El 
Asri and Hamadene [19] and El Asri [T6] extended the work of Djehiche et al [13] in the finite 



[Il[2liai6l[T5l[l9l[33lEn]. 
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horizon and infinite horizon case but allowing general jumps. In all these works the existence of 
the value functions of optimal impulse control problem and uniqueness of viscosity solution are 
obtained assuming that the switching problems without risk of default. 

Amongst the papers which consider the same problem as ours, and in the framework of 
viscosity solutions approach, the most elaborated works are certainly the ones by Djehiche and 
Hamadene |12jJ, on the one hand, and by Arnarson et al. [IJ, on the other hand. In [12j, the 
authors show existence of a solution for fll.ip . Nevertheless the paper suffers from three facts: 
(i) the switching problem have only two modes ; [ii) the switching cost functions Qij should not 
depend on x ; {Hi) the problem of switching is in finite horizon. The first issue of [12J has been 
treated by Arnarson et al. [1] since in their paper the authors show existence of the solution for 
(11.11) in the case when the growth of the functions ipi is of arbitrary polynomial type. The second 
issue of has been treated by El Asri and Hamadene [19] , since in their paper the authors show 
existence and uniqueness of the solution for fll.ip in the case when Fi = — oo. The third issu [12], 
i.e. considering the case of switching problem in infinite horizon with risk of default, was right 
now, according to our knowledge, an open problem. Note that in [Ij, the question of uniqueness 
is addressed but in the general case still remains open. Therefore the main objective of our 
work, and this is the novelty of the paper, is to show existence and uniqueness of a solution in 
viscosity sense for the systems (11. ip and (11.20 when the functions ipi, 9ij and Fi are continuous 
depending also on x and satisfy an arbitrary polynomial growth condition. We show also that 
the solution is unique in the class of continuous functions with polynomial growth. 

This paper is organized as follows: 

In Section 2, we formulate the problem and we give the related definitions. In Section 3, 
we introduce the optimal switching problem in finite horizon and infinite horizon under con- 
sideration and give its probabilistic Verification Theorem. It is expressed by means of a Snell 
envelope of processes. Then we introduce the approximating scheme which enables to construct 
a solution for the Verification Theorem. Moreover we give some properties of that solution. 
Section 4, is devoted to the connection between the optimal switching problem in finite horizon, 
the Verification Theorem and the system of PDIs (II. ip . This connection is made through back- 
ward stochastic differential equations with one reflecting obstacle in the case when randomness 
comes from a solution of a standard stochastic differential equation. We provide existence and 
uniqueness of viscosity solution of (11.10 in the class of continuous functions which satisfy a poly- 
nomial growth condition. Section 5, we show that the solution of (II. 2p is unique in the class of 
continuous functions which satisfy a polynomial growth condition. □ 
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2 Assumptions and formulation of the problem 

In finite liorizon 

Throughout this paper T (resp. A;) is a fixed real (resp. integer) positive constant. Let us 
now consider the foUowings: 

HI: b : [0,r] x IR'' IR'' and a : [0,r] x IR'' IR''^'^ are two continuous functions for 
which there exists a constant C > such that for any t G [0, T] and x, x' G 

\b{t,x)\ + \a{t,x)\ <C{l + \x\) and \a{t,x) - a{t,x')\ + \b{t,x) -b{t,x')\ < C\x - x'\ (2.1) 

H2: for i,j G T = {l,...,m}, Qij : [0,T] x ^ TR, Fi : [0,T] x M'' ^ and : 
[0,T] X iR'^ — )■ M are continuous functions and of polynomial growth, i.e., there exist some 
positive constants C and fi such that for each i,j G X: 

|^,(t,x)| + |Fi(t,a;)| + |^,,(t,x)| <C(l + |xr), V(t, x) G [0,T] xiR^ (2.2) 

H3: Moreover we assume that there exists a constant a > such that for any {t, x) G 
[0,T] X iR^ 

min{5fij(f,a;),z,j G X, i 7^ j} > a. (2.3) 

This condition means that switching from one mode to another one is not free and costs at least 
a > 0. 

We now consider the following system of m variational inequalities with inter-connected 
obstacles: V 2 G X 

mm{vi{t, x) — I max {— gij{t, x) +Vj{t,x)) V (— Fj(t,x)) 

\jei-' 

-dtVi{t, x) - Avi{t, x) - ipiit, x)} = 0, 

^ Vi{T,x) = 0, 

where X^' := X — {i} and A is the following infinitesimal generator: 

^ = i E(--)«('.-)ai-+EW*.-)|-; (2.5) 

i,j=l,k ■' i=l,k 

hereafter the superscript (*) stands for the transpose, Tr is the trace operator and finally < x,y > 
is the inner product of x, y G IR''. 

The first main objective of this paper is to focus on the uniqueness of the solution in viscosity 
sense of (12. 4p . To proceed we will precise the notion of a viscosity solution of the system (12. 4p . 
It will be done in terms of subjets and superjets. 



Definition 1 Let v G C((0,T) x iR^), {t,x) an element o/(0,T) x IR^ and finally Sk the set of 
k X k symmetric matrices. We denote by J'^'^v(t,x) (resp. J'^'^v{t, x)), the superjets (resp. the 
subjets) of V at {t, x), the set of triples {p, q, X) E M x x Sk such that: 

v{s, y) < v{t, x) + p{s -t) + {q,y - x) + l{X{y - x),y - x) + o{\s - t\ + \y - x\^) 
{resp. v{s, y) > v{t, x) + p{s — t) + {q,y — x) + ^{X{y — x),y — x) + o(|s — t\ + \y — x\'^)).n 

Note that if ip — v has a local maximum (resp. minimum) at {t, x), then we obviously have: 

{Dtip{t,x),D^ip{t,x),Dl^ip{t,x)) e J^'-v{t,x) (resp. J'^'+v{t,x)).D 

We now give the definition of a viscosity solution for the system of PDE equations with risk 
of default in finite horizon fl2.4p . 

Definition 2 Let {vi, ...,Vm) be a m-uplet of continuous functions defined on [0,T] x , M- 
valued and such that {vi, ...,Vm)(T,x) = for any x G IR'' . The m-uplet {vi,...,Vm) is called 
a viscosity supersolution (resp. subsolution) of Iji2.4\ ) if for any i G I, {t, x) G (0, T) X m}" and 
{p,q,X) G J'^''Vi(t,x) (resp. J^'-^ Vi{t , x) ) , 

mm < Vi{t, x) — I max {— gij{t, x) + Vj{t, x)) V {—Fi{t, x)) 

-p - lTr[a*Xa] - {b, q) - %l)i{t, x)} > {resp. < 0). 
It is called a viscosity solution it is both a viscosity subsolution and supersolution .□ 

As pointed out previously we will show that system f l2.4l) has a unique solution in viscosity sense. 
A particular case of this system is the deterministic version of the optimal m-states switching 
problem in finite horizon with risk of default which is well documented e.g. in [H |T2] and which 
we will describe in the next section. 
In infinite horizon 

Let us now consider the followings assumption: 

H4: h : ^ IR^ and a : — )■ IR^^'^ are two continuous functions for which there exists a 
constant C > such that for any x, x' G IR^ 

|fo(x)| + |a(a;)| < C(l + and \a{x) - a{x')\ + \b{x) - b{x')\ < C\x - x'\ (2.6) 

H5: for i,j G X = {1, ...,m}, Qij : IR^ — )■ iR is a continuous function. Moreover we assume 
that there exists a constant a > such that for any x G M'^, 

-<giA^)<a, ^1,3 ex, i^j. (2.7) 
a 
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H6: for i & X ipi : — > IR and Fj : — > IR^ are continuous function of polynomial 
growth, i.e., there exist a constant C and ^ such that for each i eX: 

\iPi{x)\ + < C(l + Vx e IR\ (2.8) 

We now consider the following system of m variational inequalities with inter-connected 
obstacles: V 2 G X 



mm 



Vi{x) - imax{-gij{x) + Vj{x)) V (-Fj(x)) j ,rvi{x) - Avi{x) - tpiix) | = 0' (2-9) 



where X * := X—{i}, r is a positive discount factor and A is the following infinitesimal generator: 

ij=l,k ■' i=l,k 

The second main objective of this paper is to focus on the uniqueness of the solution in 
viscosity sense of (12. 9p . We now give the definition of a viscosity solution of the elliptic system 
with inter-connected obstacles (12. 9p . 

Definition 3 Let {vi,...,Vm) he a m-uplet of continuous real-valued functions defined on . 
The m-uplet (vi, ...,fm,) called a viscosity supersolution (resp. subsolution) of 1^2. 4\ ) if for any 
i eX, X e IR^ and {q,X) e ,P'^Vi(t,x) (resp. J'^'^Vi{x)), 

min |vj(x) - |^max(-fi^y(x) + Vj{x)) V (-Fi(x))^ , 

rvi{x) - \Tr[a*Xa] - {b, q) - V^i(x)} > {resp. < 0). 
It is called a viscosity solution if it is both a viscosity subsolution and supersolution . □ 

As pointed out previously we will show that system (12. 9p has a unique solution in viscosity 
sense. This system is the deterministic version of the optimal m-states switching problem in 
infinite horizon with default risk which is well documented in [T3| [T7] and which we will 
describe briefly in the next section. 

3 The optimal m-states switching problem 
3.1 In finite horizon with risk of default 

Let (i7, J-", P) be a fixed probability space on which is defined a standard d-dimensional Brownian 
motion B = {Bt)o<t<T whose natural filtration is (J-'° := (j{Bs,s < t})o<t<T- Let F = {J^t)o<t<T 
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be the completed filtration of (J^°)o<t<T with the P-nuU sets of hence (J^t)o<t<r satisfies the 
usual conditions, i.e., it is right continuous and complete. Furthermore, let: 

- be the cr-algebra on [0, T] x Q of F-progressively measurable sets; 

- JkA?'^ be the set of T'-measurable and TR*^- valued processes w = {wt)t<T such that E[j^ \ tVs\^ds\ < 
oo and <S^ be the set of P-measurable, continuous processes w = {wt)t<T such that E[s\iY)^^rp p] < 
oo; 

- for any stopping time r G [0,T], %. denotes the set of all stopping times 9 such that 
T<e<T. 

The problem of multiple switching can be described through an example as follows. Assume 
we have a plant which produces a commodity, e.g. a power station which produces electricity. 
The production activity have m modes, or "definitely closed/defaulting" indicated by f. A 
management strategy of the plant consists, on the one hand, of the choice of a sequence of 
nondecreasing stopping times (r„)„>i {i.e.Tn < Tn+i and tq = 0) and the stopping time 7 where 
the manager decides to switch the activity from its current mode to another one or definitely stop 
the production. On the other hand, it consists of the choice of the mode a r.v. J>„-measurable 
with values in I, to which the production is switched at t„ from its current mode. Therefore the 
admissible management strategies of the plant are the pairs (5,^) :— ((t„)„>i, 7, (^„)„>i) and 
the set of these strategies is denoted by T>. 

Let now X :— {Xt)o<t<T be an T'-measurable, M'^-valued continuous stochastic process which 
stands for the market price of k factors which determine the market price of the commodity. On 
the other hand, assuming that the production activity is in mode 1 at the initial time t — 0, let 
{ut)t<T denote the indicator of the production activity's mode at time t G [0,T] : 

Ut = l[0,ri] 

n>l 

Then for any t <T, the state of the whole economic system related to the project at time t is 
represented by the vector : 

(t, Xt, ut), if Tn < t < r„+i; 
(7, X^), if in mode f . 

Finally, let 'i/ji{t,Xt) be the instantaneous profit when the system is in state {t,Xt,i), for 
i,j & I i ^ j, let gij{t, X^) denote the switching cost of the production at time t from current 
mode i to another mode j and let Fi{'y,Xj) denote the cost of default (definitely stop the 
production) at time 7, when in mode i and denote -Fi(7, X^) — F{'y, X^, u^) when — i. Then 
if the plant is run under the strategy (S,^) — ((t„)„>i, 7, (^„)„>i) the expected total profit is 
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given by: 

•^0 n>l 

Therefore the problem we are interested in is to find an optimal strategy, i.e, a strategy {S*^u*) — 
((<)n>i,7*), (O) such that J{5*,C) > mo for any {5,0 e V. 

Note that in order that the quantity J{S,^) makes sense we assume throughout this paper 
that for any i,j el the processes {Fi{t, Xt))t<T , {gij{t, Xt))t<T (resp. {ipiit, Xt))t<T) belong to 

(resp. A^^'^). On the other hand there is a bijective correspondence between the pairs (5, 
and the pairs {S,u). Therefore throughout this paper one refers indifferently to {5,^) or {5,u). 

The verification Theorem for the m-states optimal switching with risk of default problem is 
the following: 

Theorem 1 Assume that there exist m processes {Y^ :— (^*)o<t<r,^ = 1, ■■■,m) of such that: 
yt < T, Yi = ess sup^^tE[j;iji{s,Xs)ds + max(-^i,(r, + Yj) V Fi{T, Xr)l[r<T]\J^t], 

(3.3) 

Then: 

(i) Yq = sup J {5, u). 

{ii) Define the sequence of ¥ -stopping times 6 = {tI,j), i = 1,2, ...,m, < t < T. as follows : 

ri = al^^^T, 

where: 

• The first time the activity defaults while in mode i is given by 

a\ :=inf{s >t, Y^' = Fi{s, Xs)} AT, i = l,...,m. 

• The first time the activity is switched from mode i to any of the other modes j ^ i is 
given by 

a\ := inf{s > t, Yl = max(-5,,(s, X^) + F/)} A T 
Finally, let j = sup r^*. Then, the strategy {6= {{Tt)t>o,7),u*) is optimal. □ 

0<t<T 
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Proof. The arguments of proof are standard, based on the properties the Snell envelope and is 
proved in [13], Theorem 1. □ 

The issue of existence of the processes Y^,...,Y^ which satisfy fl3.10p is also addressed in 
[13] . Also for n > let us define the processes [Y^'^, ...,Y"^'^) recursively as follows: for i eX 
we set, 

y/'° = ess snVr>tE[j^ Ms,X,)ds + Fiir, Xr)l[r<T]\J'tl 0<t<T, (3.4) 

and for n > 1, 

y/'" = ess sup,>,^[ / ilJiis,X,)ds + max{-gikir,Xr) + Y^'^''-')VF,{r,Xr)l[r<T]\J't], 0<t<T. 

(3.5) 

Then the sequence of processes ((y^'", y"^'"'))„>o have the following properties: 

Proposition 1 (^^, Pro. 3 and Th.2) 

{%) for any z G X and n > 0, the processes y^'", ...,y™'"' are well-posed, continuous and belong 
to S^, and verify 

Wt < T, Yr < rr^' < E[ Tjmax \^i{s, X,)\}ds\J^t]; (3.6) 

Jt *=i'™ 

(ii) there exist m processes Y^ , y™ of such that for any i G X; 
(a) Vt < T, Yi = lim / y/'" and 



E[sup |y;'" - y;r] ^ O as n^+oo 

s<T 



(b) Vt < T, 

Yi = ess sup,^tE[ [ ij^is, Xs)ds+maxi-g,kiT, Xr) + Y^^)y F^ir, Xr)l[r<T]\J't] (3.7) 
Jt '^'^^ ' 

i.e. y^,...,y™ satisfy the Verification TheoremlMo 



Remark 1 Note that the characterization IjjS. 7| ) implies that the processes Y^, Y^ of which 
satisfy the Verification Theorem are unique. 

3.2 In infinite horizon with risk of default 

Let {Vt, J-", P) be a fixed probability space on which is defined a standard d- dimensional Brownian 
motion B = {Bt)t>o whose natural filtration is (7-"° := a{Bs,s < t})t>o. Let F = {J-'t)t>o be 
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the completed filtration of {J^)t>o with the P-nuU sets of J^, hence {Tt)t>o satisfies the usual 
conditions, i.e., it is right continuous and complete. Furthermore, let: 

- be the cr-algebra on [0, +00) x Q of F-progressively measurable sets; 

- M'^''^ be the set of T'-measurable and Ft'^-valued processes w — {wt)t>o such that 
E[Io°° I'^sl'^ds] < 00 and be the set of P-measurable, continuous processes w — {wt)t>o such 

that i?[supj>o l""^*!^] < '^1 

- for any stopping time r G M^, % denotes the set of all stopping times 6 such that t < 6; 

- for any stopping time r, J> is the cr-algebra on fl which contains the sets A of J-' such that 
An{T <t} e J^t for every t > O.D 

A decision (strategy) of the problem of multiple switching, on the one hand, consists of the 
choice of a sequence of nondecreasing stopping times (T„)n>i {i-e.Tn < r^+i and Tq = 0) and 
the stopping time 7 where the manager decides to switch the activity from its current mode to 
another one or definitely stop the production. On the other hand, it consists of the choice of 
the mode a r.v. J>^-measurable with values in I, to which the production is switched at 
Tn from its current mode. Therefore the admissible management strategies of the plant are the 
pairs {6,$,) :— ((t„)„>i,7, (Cn)n>i) and the set of these strategies is denoted by V. 

Let now X :— {Xt)t>o be an T'-measurable, l?*^-valued continuous stochastic process which 
stands for the market price of k factors which determine the market price of the commodity. On 
the other hand, assuming that the production activity is in mode 1 at the initial time t — 0, let 
('Wt)t>o denote the indicator of the production activity's mode at time t e : 

ut = l[0,ri](i) + ^-l(-n,r„+i](i)- (3.8) 

n>l 

Then for any i > 0, the state of the whole economic system related to the project at time t is 
given by the vector: 

{t,Xt,Ut), \lTn<t<Tn+l] 

(7, X^), if in mode f • 
Finally, let ijJiiXt) be the instantaneous profit when the system is in state {t,Xt,i), for 
i,j G X i ^ j, let gij{Xt) denote the switching cost of the production at time t from the 
current mode i to another mode j and let -Fj(7,X^) denote the cost of default (definitely stop 
the production) at time 7, when in mode i and denote Fi{X^,) = F{X^,u^) when = i. When 
the plant is run under the strategy (5, ^) = ((r„)„>i, 7, (^n)n>i) the expected total profit is given 
by: 

n>l 
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Then the problem we are interested in is to find an optimal strategy, i.e, a strategy [S*,^*) such 
that J{6*,C) > J{S,0 for any (5,0 e V. 

Note that in order that the quantity J{6, ^) makes sense we assume throughout this paper 
that for any i E X the processes {e~'^^tl)i{Xt))t>Q and (e~''*Fj(Xt))t>o belong to M."^'^ and S"^ 
respectively. 

The Verification Theorem for the m-states optimal switching with risk of default problem in 
infinite horizon is the following: 

Theorem 2 . Assume that there exist m processes (F* := {Y^)t>Q,i = 1, ...,m) of such that: 

yt > 0, e-^*y/ = ess sup,^tE[j^ e-'-'ij,{Xs)ds + e"^" max{-gij{Xr) + Yj) V e-^"F,(X,)| J^^], 
lim (e-'^*y;) = 0. 

(3.10) 

Then: 

(i) Iq*^ = sup J {6, u). 

(a) Define the sequence of F-stopping times 5 = (r^*,7), i = 1,2, ...,m, t >0. as follows : 
where: 

• The first time the activity defaults while in mode i is given by 

ai := inf{s > t, = F,(X,)}, i = 1, ...,m. 

• The first time the activity is switched from mode i to any of the other modes j ^ i is 
given by 

al := inf{s > t, Y^ = max(-^,,.(s, X,) + F/)} 

Finally, let 7 = supr/. Then, the strategy (5 = ((Tt)t>o,7) is optimal. □ 

t>o 

Proof. The arguments of proof are standard, based on the properties the Snell envelope and is 
proved in [TB], Theorem. 1. □ 

The issue of existence of the processes Y^,...,Y'^ which satisfy (13.101) is also addressed in 
[13] . For n > let us define the processes (F"'^, ...,Y^'"^) recursively as follows: for z G X we 
set, 

^-rtyOA ^ SUp,>,F[^ e-"V^i(X,)ds + 6-'^ Fi{Xr)\Tt], t > 0, (3.11) 
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and for n > 1, 

g-rtyn,. ^ sup,>,^[ [ e-'-'iJi{Xs)ds + e-'^ max(-(/,fc(X.) + y;^-i''=) Ve-^^F,(X,)|J-i], t > 0. 

(3.12) 

Then the sequence of processes {{Y^'^ , ^"''™))n>o have the following properties: 

Proposition 2 f/I^, Pro. 3 and Th.2) 

{i) for any i eX and n > 0, the processes Y"^'^ , ...,y"'™ are well-posed, continuous and belong 
to S^, and verify 

r+oo 

Vt > 0, e-''*F"'' < e-''*y"+''' < £;[ / e-'-'{imix\i;,{X,)\}ds\J^t\] (3.13) 

Jt *=i'™ 

(ii) there exist m processes , Y"^ of such that for any i G X: 

(a) Vt > 0, Y; = lim / y"'^ 

n— >oo 

(b) Vt > 0, 

e~^*y/ = ess sup,>,E[ [ e-"^i(X,)ds + e"^^ max (-(?,, (X,) + y,'^) V e-^^F,(X,)| J-^] 

(3.14) 

i.e. Y^,...,Y"^ satisfy the Verification TheoremlM^ 



Remark 2 The characterization Iji3.14 ) implies that the processes Y^, ...,y™ of which satisfy 



the Verification Theorem are unique.D 

4 Uniqueness of the viscosity solution in finite horizon 

Let {t,x) G [0,T] X and let (X*^)s<r be the solution of the following standard SDE: 

t^Xf = b{s, Xf)ds + a(s, Xf)dBs for t < s < T and Xf = x for s < t (4.1) 

where the functions b and a are the ones of (12.11) . These properties of a and b imply in particular 
that the process (X*^)o<s<t solution of the standard SDE (14.11) exists and is unique, for any 
t G [0, T] and x G iR^ 

The operator A that is appearing in (12. 5p is the infinitesimal generator associated with X*'^. 
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Theorem 3 There are deterministic functions v^, ...,v"^ : [0, T] x — > M such that: 

y{t,x) e [o,r] X iR^vs G [t,r],y;'*^ = v'{s,xl^), i = i,...,m. 

Moreover the functions {v^, ...,v'^) : [0,T] x — )■ M are continuous, solution in viscosity sense 
of the system of variational inequalities with inter- connected obstacles ^^K^and of polynomial 
growth.D 

Proof: The continuity of functions v^,...,v"^ foUows from the dynamic programming principle 
and is proved in [19]. □ 

Now we give an equivalent of quasi- variational inequality (12 ■4p . In this section, we consider 
the new function Fj given by the classical change of variable rj(t,x) = exp{t)vi{t, x), for any 
t G [0,r] and x G iR^. Of course, the function Fj is continuous and of polynomial growth with 
respect to its arguments. 
A second property is given by the 

Proposition 3 Vi is a viscosity solution of Iji2.4\ ) if and only if Fj is a viscosity solution to the 
following quasi-variational inequality in [0,T[xlR^ , 

mm{Ti{t,x) - ( max(-e*5fy(t,x) + Tjit,x)) V (-e*Fi(t, x)) ) , 
Vie!-' / 

Ti{t,x)-dtTi{t,x)-ATi{t,x)-e%{t,x)} = 0, (^'2) 

^ Ti{T,x) = e^v,{T,x) = 0.n 

We are going now to address the question of uniqueness of the viscosity solution of the system 
(12.41) ■ We have the following: 

Theorem 4 The solution in viscosity sense of the system of variational inequalities with inter- 
connected obstacles ^2.4^ is unique in the space of continuous functions on [0,T] x which 
satisfy a polynomial growth condition, i.e., in the space 

C := {(f : [0, T] X IR^ — ?■ IR, continuous and for any 

{t, x), \(p{t, x)\ < C(l + l^l^) for some constants C and /i}. 

Proof. We will show by contradiction that if ui, ...,Um and wi, ...,Wm are a subsolution and a 
supersolution respectively for (14. 2 p then for any i = 1, ...,m, Ui < Wi. Therefore if we have two 
solutions of ( 14. 2 p then they are obviously equal. Actually for some R > suppose there exists 
{t,x,J) G (0,T) X BrxI {Br := {x G R''; \x\ < R}) such that: 

max{ui(t, x) — Wi{t, x)) = uj(i,x) — wj(t,x) = ri > 0. (4.3) 
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Let us take 6, A and /3 G (0,1] small enough, so that the following holds: 



/3T<| 

-\w-{lx)<l (4.4) 

t ^ 4- 

Here 7 is the growth exponent of the functions which w.l.o.g we assume integer and > 2. Then, 
for a small e > 0, let us define: 

^%x,y) = u,{t,x)-{l~\)w,{t,y)-^Jx-y\'''-ei\x-x\'''-'' + \y-x\'^^^^ (4.5) 

By the growth assumption on Ui and Wi, there exists a {to, xo,yo,io) G (0,r] x Br x Br x X, 
such that: 

^'^'{to,Xo,yo) = max <l>*(t,x,y). 

On the other hand, from 2^1° {to, Xo, yo) > '^'^'(^o, Xo, xo) + $e°(to, ^/o, yo), we have 

- ?/oP'^ < {uu,{to,Xo) - u,,,{to,yo)) + (1 - X){wi^{to,xo) - Wi„(to, l/o)), (4.6) 

and consequently ^|xo — is bounded, and as e — )■ 0, l^o — yo\ — ?■ 0. Since and Wi^ are 

uniformly continuous on [0,T] x S/j, then ^|xo — — )■ as e — )■ 0. 

Since 

(1 - \)uj(t,x) - wj(t,x) - 7 < $^"(to,a:o,l/o) < (1 - A)^Zio(^o, a:o) - Wi^{t^,yo) - ^, 

t to 

it follow as A —7- and the continuity of u and w that, up to a subsequence, 

{to,Xo,yo,io) ^ {t,x,x,i). (4.7) 
Next let us show that to < T. Actually if to = ^ then, 

^,{t,x,x)<^i'{T,xo,yo), 

and, 

uj{t,x) - (1 - X)wj{t,x) - ^ < -/3(T - If - ^, 
since Uif^{T,xo) = Wj„(T, ?/o) = 0. Then thanks to (14.31) we have, 

A 



V <-Xwj{t,x)+PT + ^ 



V < zV- 

which yields a contradiction and we have to ^ (0)^)- We now claim that: 

Ui^Xto,xo) - ( max {-e'^''gioj{to,Xo) +Uj{to,xo)} V {-e^" Fi^{to, xo)) ] > 0. (4. 
\jex-^o 
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Indeed if 

Uio{to,xo)- I max {-e^° gi,,j{to, xo) + Uj{to,xo)}\/ {-e^"Fi^,{to,xo))] < 0. 

Casel : max {-e^''gioj{to,Xo) + Uj{to,Xo)} V {-e*° Fi^{to, Xq)) = -e*° Fi^{to,Xo). 
Then 

Uio{to,xo) < -e*"Fio(to,a;o). 
From the supersolution property of Wig(to, i/q), we have 

Wi^{to,yo) > I max {-e^°gi^j{to,yo) + Wj{to,yo)} y {-e^°Fi^{to,yo)) 
then 

Wi^{to,yo) > -e*"Fio(to,yo)- 

It follows that: 

Uig{to,xo) -Wi^{to,yo) < -e^^iFi^ih.xo) - Fi^{to,yo)). 

But this contradicts the definition of (14.31) . since F, u, w is uniformly continuous on [0,T] x Bji 
and sending A — )■ and the claim (14. 8 p holds. 

Case2 : max {-e^ogi^j{to,xo)+Uj{to,Xo)}y{-e^''Fig{to,xo)) = max {-e^°gigj{to,xo)+Uj{to,xo)}. 
Then there exists k G X^*" such that: 

UiQ{to,xo) < -e^°giok{to,xo) + Uk(to,xo). 
From the supersolution property of Wig{to, yo), we have 

Wioito, yo) > ( max {-e^^giojih, yo) + Wj{to, yo)} V {-e^'> Fi^{to, yo)) ) , 

then 

Wi^,{to, yo) > -e^°giok{to, yo) + Wfc(to, Vo)- 

It follows that: 

Ui^{to,xo)-il-X)wio{to,yo)-{uk{to,xo)-{l-X)wkito,yo)) < {l-X)e^''giok{to,yo)-e*°9iok{to,xo). 

Now since gtj > a > 0, for every i ^ j, and taking into account of (14.51) to obtain: 

¥°{to,xo,yo) - $^(to,a;o,2/o) < -aAe*« + giok{to,yo) - e*°giok{to,Xo). 

But this contradicts the definition of iq, since gi^k is uniformly continuous on [0,T] x and 
the claim ( Oil holds. 
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Next let us denote 



Mt, x,y) = j-\x- y\'^ + e{\x- ^ 1^ _ -|2,+2) ^ -tY + ^. (4.9) 

Then we have: 

Dt^,{t,x,y)=2(]{t-t)-^, 

D^ip^it, X, y) = -y)\x- y\^''-^ + 9(2^ + 2)(x - x)\x - x\'\ 
Dyif^it, X, y) = -l{x- y) \x - yp-'' + ^(27 + 2)(y - x)!?/ - x\^\ 

ai{x,y) -ai{x,y)\ ^ (a2{x) \ (4.10) 
-ai{x,y) ai{x,y) J y a2{y) J 

with ai{x, y) = ^\x — yp'''"^/ + 7(27 — 2)(x — ?/)(x — ?/)*|x — y\'^^~'^ and 
a2(x) = ^(27 + 2)|x - xpTJ + 2^7(27 + 2)(x - x){x - xy\x - 

Taking into account f l4.8p then applying the result by Crandall et al. (Theorem 8.3, [7]) to the 
function 

Ui^,{t,x) - (1 - \)wi,,{t,y) - ip,{t,x,y) 
at the point (to, Xq, yo), for any ei > 0, we can find c,d E IR and X, 1" G 5'^, such that: 



B{t,x,y) = Dl ip,{t,x,y) 



(c, f (xo - 2/o)ko - yoP^-' + ^(27 + 2){xo-x)\xo - x\^\X) G j2'+(«,„(to, xq)), 
(-rf, ^(xo - yo)\xo - yo^'^' - 0{2j + 2){yo -x)\yo -x\'\Y) G J2'-((l - X)w,,{to,yo)), 
c + d = Dtifeih, xo, yo) = 2;5(to -t) - ^ and finally 

\ 

j^ + \\B{to,Xo,yo)\\)I < I ^ ^1 < B{to,xo,yo) + eiB{to,Xo,yof. 

(4.11) 

Taking now into account (14. Sp . and the definition of viscosity solution, we get: 

-c + Ui,Xto,xo) - ^Tr[a*{to,xo)Xa{to,xo)] - {^{xq - yo)\xo -?/oP^"^+ 

6'(27 + 2)(xo — x)\xo — xp^, 6(to, a;o)) — e'^°i^io(io, Xq) < and 

d+{l- X)wi,{to, yo) - ^Tr[a*{to, yo)Ycr{to, yo)] - (^(xq - yo) \xo - ^/oP^"'- 

9i2^ + 2)iyo-x)\yo-x\'\b{to,yo)) - (1 - A)e*°^,o(to, 2/0) > 
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which imphes that: 

-c - d + Uio{to,xo) - (1 - X)wi^{to,yo) < lTr[a*{to, xo)Xa{to, xq) - a*{to,yo)Ycr{to,yo)] 

+ (7(3^0 - yo)\xo - ^/oP'^"^&(^o,a;o) - &(^o,?/o)) 
+ (^(27 + 2){xo - x)\xo - x\^\ b{to, xo)) 
+ {e{2j + 2)iyo - x)\yo - x\^^, h{U, yo)) 
+e''>iJi,{to,xo) - (1 - A)e*«^/;i,(fo,2/o). 

(4.12) 

But from fl4.10p there exist two constants C and Ci such that: 

\\aiixo,yo)\\ < C\xo - yol'""' and (||a2(xo)|| V IMoW < 0,6. 

As 

„ N 1 I «i(a;o,2/o) -oi(a;o,?/o) \ / a2(a;o) 

B = B{to,xo,yo) = - \ ^ + 

e \-ai{xo,yo) ai(xo,?/o) / \ 02(2/0) 



then 

< 

e " \ -/ / 

It follows that: 



B<-\xo-yo\''-'\ ^ ^ \+C,9I. 



B + t,B^ < C{-^\x, - yoP^-' + |l^o - 2/0^^-') y ^ I + C^i^/ (4.13) 

where C and Ci which hereafter may change from line to line. Choosing now ei = e, yields the 
relation 

B + e,B^ < ^(|xo - yor^-' + 1^0 - y,\''-') ^ ^ ~l I + C^ei. (4.14) 
Now, from ([211]), fHTT]) and (gUD we get: 

^rr[a*(to,xo)Xa(to,xo)-a*(to,2/o)l^a(to,2/o)] < -(|a;o-?yo|'^+|a;o-2/o|'^-')+Cie(l+|xop+h- 



Next 



and finally, 



e e 



(^(27 + 2){xo - x)\xo - x\^\ b{to, Xo)) + {e{2-f + 2){yo - x)\yo - xp^ b{to, yo)) 

< 90(1 + \xo\\xo - xp^+i + bolbo - 5^1'^+'; 
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So that by plugging into f l4.12p and note that A > we obtain: 

-2^(to - t) + I + Ui,{to, xo) - (1 - X)w^,{to, yo) < f (|xo - ?/oP^ + \xo - 2/o|'^-') + 

Ci9{l + |xop + boH + ^|xo - yo\'^ + 0C{1 + \xo\\xo - x\^^+' + \yo\\yo - x\'^^')+ 
e^°i:io{to,xo) - (1 - A)e*«V^io(to,?/o)- 
By sending e— J-O, A— J-O, ^->0 and taking into account of the continuity of i/ji^ and 7 > 2, we 
obtain i] < which is a contradiction. The proof of Theorem H] is now complete. □ 

As a by-product we have the following corollary: 

Corollary 1 Let {v^, ...,v"^) be a viscosity solution of i2.4\) which satisfies a polynomial growth 
condition then for i = 1, ...,m and (t, x) G [0, T] x IR^ , 

v\t,x) = sup E[ ri;uSs,Xf)ds -J^gu^^ (r„, X*^)l[,„<^] - F(7, i/^)l[^<T]].n 

5 Uniqueness of the viscosity solution in infinite horizon 

Let X G M'^ and let be the solution of the following standard SDE: 

dX^ = b{X^)dt + (r{X^)dBt, X^ = X (5.1) 

where the functions b and a are the ones of H4. These properties of a and b imply in particular 
that X^ solution of the standard SDE (15.11) exists and is unique in ]R^. The operator A defined 
in (I2.10p is the infinitesimal generator associated with X^. 

Theorem 5 There are deterministic functions v^, ...,v"^ : M!^ — ?■ M such that: 

yx e R^Y^''' = v\x), 2 = l,...,m. 

Moreover the functions {v^, ...,t'™') : iR^ — > M are continuous, solution in viscosity sense of the 
system of variational inequalities with inter- connected obstacles ^2. Pj) and of polynomial growth. 

Proof: The continuity of functions w^, follows from the dynamic programming principle 

and is proved in [16]. □ 

We are going now to address the question of uniqueness of the viscosity solution of the system 
(12.91) . We have the following: 

Theorem 6 The solution in viscosity sense of the system of variational inequalities with inter- 
connected obstacles l{2.9\) is unique in the space of continuous functions on which satisfy a 
polynomial growth condition, i.e., in the space 

C := {ip : — )■ ]R, continuous and for any 

X, |v^(2;)| < (^(1 + \x\^) for some constants C and /i}. 
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Proof : We will show by contradiction that if ui, ...,Um and Wi, ...,'Wm are a subsolution and a 
supersolution respectively for fl2.9p then for any i = 1, ...,m, Ui < Wi. Therefore if we have two 
solutions of (12. 9 p then they are obviously equal. Actually for some R > suppose there exists 
{xq, io) e BrX Z {Br ■= {x e M''; \x\ < R}) such that: 

max{ui{x) - Wi{x)) = Uig{xo) - Wig{xo) = r] > 0. (5.2) 

Then, for a small e > 0, and ^, A G (0, 1) small enough, let us define: 

$:(a;, y) = u^ix) - (1 - X)w,iy) - ^1^ - vl'^ - 0{\x - xoP^+' + \y - xo\'''+'). (5.3) 

By the polynomial growth assumption on Ui and Wi, there exists a {x^^ye^ie) £ Br x Br x X, 
such that: 

$Hx„?/,) = m&x^\{x,y). 
On the other hand, from 2$*'(xe, Ty^) > <l>*'(xe,Xe) + $*'(ye,ye), we have 

Y^e - Vel'^^ < {UiS^e) " UiSVe)) + (1 " A)(u)i^(x,) - Wi^y,)) 

< ^ \Ui{Xe) - Ui{ye) \ + (1 - A) X] Iwiixe) - Wi{ye)\ ^^'^^ 

and consequently ^\x^ — y^\'^'^ is bounded, and as e — > 0, jx^ — ye| — > 0. Since Ui and Wi are 
uniformly continuous on Br, then ^jx^ — yel"^^ — )■ as e — )• 0. 
Since 

Mio(xo) - (1 - X)wi^ixo) < ^';{x„y,) < Ui^{x,) - (1 - X)wi^{ye), 
it follow as A —7- and the continuity of Ui and Wi that, up to a subsequence, 

(x,,y,,z,) (xo,Xo,zo)- (5.5) 

We now claim that: 

UiXxe) - i max {-gi^j{x^) +Uj{x^)} V {-Fi^{x^)) ) > 0. (5.6) 

Indeed if 

Ui^{x^) - ( max {-Qi^jix^) + Uj{x^)} V (-Fj^(x,)) ) < 0. 

Casel : ( max {-gi,j{x,) + Uj{x^)} V {-Fi^{xJ) ) = -Fi^{x,). 
Then 

Wi,(a;e) < -Fi^{xe). 
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From the supersolution property of WiXVe)-, we have 



WiXVe) > max {-gijive) + Wj{y,)} V {-FiXVe)) 
then 

It follows that: 

But this contradicts the definition of f l5.2p . since F, w is uniformly continuous on 5^ and 
sending A — )■ and the claim (15.61) holds. 

Case2 : I max {-gij{x^) + Uj{x^)} V (-Fj^(x,)) ) = max {-gij{x^) + Uj{x^)}. 
then there exists k G X~*" such that: 

From the supersolution property of Wi^y^), we have 



WiXVe) > max {-gijiVe) + Wj{ye)} V {-FiXye)) 
then 

It follows that: 

Ui^Xe) - (1 - A)Wi^(?/,) - (Mfc(Xe) - (1 - A)tyfc(2/e)) < (1 " >^)giAye) " fi'i.fcla^J- 

Now since gij > a > 0, for every i 7^ j, and taking into account of fl5.3l) to obtain: 

¥;{x,,y,) -<l>^(x,,y,) < -aX + gi^y,) - giA^e) 

But this contradicts the definition of i^, since gi^k is uniformly continuous on Br and the claim 
(EH) holds. 

Next let us denote 

ip^x, y) = l|x - + e{\x - xoP^+' + \y- xoP^+'). (5.7) 
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Then we have: 



Dyip,{t,x,y) 



^{x -y)\x- + e{2j + 2){x - xo)\x - xo\^^, 

-y)\x- i/|27-2 + 0^2^ + 2){y - yo)\y - yo\'\ 



< B{t,x,y) = Dlyip,{t,x,y) 



1 j aiix,y) -ai{x,y) 
-ai{x,y) ai{x,y) 



a2{x) 
(32 (y) 



(5. 



with ai(x, y) = ^\x - y\'^'^ ^/ + 7(27 -2){x - y){x - y)*\x - y\^"< ^ and 
02(2;) = ^(27 + 2)\x- XoP^/ + 2^7(27 + 2)(x - Xo){x - Xoy\x - XqP^"^ 

where I stands for the identity matrix of dimension k. Taking into account (I5.6p then applying 
the resuh by Crandall et ah (Theorem 3.2, [7j) to the function 

Ui{x) - (1 - X)Wi{y) - Lp^{x, y) 



at the point (xe, ye), for any ei > 0, we can find X, F G Sk, such that: 



(^(xe - ye)\xe - ye\'^-' + 0(27 + 2)(x, - xo)|x, - xo\^\ X) G J2'+K(x,)), 
C-{xe - ye)\xe - ye\'^-' - ^(27 + 2){ye - yo)\ye - yo\'\ Y) G - X)wdy.)), 

-(i + ||S(a;„?/,)||) °j < <B{xe,ye) + e,B{xe,ye?. 

Taking now into account (15. 6p . and the definition of viscosity solution, we get: 

ru^Xxe) - \Tr[a*{xe)Xa{xe)] - (x, - - y.p^-^ 

+^(27 + 2){x, - xq)\x, - xop^, h{xe)) - ^iX^e) < and 
r(l - X)wUye) - lTr[a*{ye)Ya{ye)] - (^(x, - - y.p^-s 

-^(27 + 2){ye - xo)\ye - Xo\^\ Kz/e)) - (1 - X)i^dy.) > 

which implies that: 

rui^xe) - r(l - X)wUye) < |Tr[a*(x,)Xa(x,) - a*iye)Ya{ye)] 
+ C-{xe - ye)\xe - y^r'-', b{xe) - b{ye)) 

+ {d{2-f + 2){xe - xo)\xe - x^\^\ b{xe)) + (^(27 + 2)(y, - xo)\ye - xqP^ 

+ ^^JiSXe)-{l-^)^^Xye)■ 

But from (15. 8p there exist two constants C and Ci such that: 

\\ai{xe,ye)\\ < C\xe - ye\^^-^ and (||a2(x,)|| V IMyM < C.O. 
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As 



then 



1 / ai{xe,ye) -ai{x^,y,)\ /a2(x,) 
B = B{x„ye} = - I + 

-ai{x^,y,) ai{x^,y,) \ a2(2/e) 



B<-\x,-y,\'''-' I ^ J +^^^0 / 



It follows that 

,27-2 , £l 



where C and Ci which hereafter may change from line to line. Choosing now ei = e, yields the 
relation 

Now, from H4, ([53]) and fl5J2|) we get: 



Next 



-(x, - yJla;. - ye?'^ ^ ^(a^e) - Ky^)) < — \x, - y, 
e e 

and finally, 

(0(27 + 2)(Xe - Xo)|Xe - Xo\^\ h{x,)) < 6*0(1 + \x,\)\x, - Xo|'^+' 
(0(27 + 2){y, - x,)\y, - x,\^\ h{y,)) < eC{l + \y,\)\y, - x,\^^^\ 
So that by plugging into flS.lOp we obtain: 

rU^X:) - r(l - \)w,Xy.) < f (l^e - + W - + C^e{l + + + 

- y.?' + ^C'(l + \xX)\x, - xo|2^+i + eC{l + ly.DI?;, - a;oP^+i + 

'IpiX^e) - (1 - A)^/'i^(?/,). 

By sending e— i-O, A— i-O, 0— J-O and taking into account of the continuity of -i/'i^, we obtain 
UiQ{xQ) — t0jo(2;o) < which is a contradiction. The proof of Theorem [6] is now complete. □ 
As a by-product we have the following result: 

Corollary 2 Let (f^, ...,v^) he a viscosity solution of Ii2.9\) which satisfies a polynomial growth 
condition. Then for i = I, ...,m and {t,x) G M'' , 

sup E[ f e-^^^uXX:)ds - $^e— ^„,JX^J - e-^^F(X^^ 



V {x} = 

(5,0 eoj -^0 



n>l 
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